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ABSTRACT 


\> 

Exact  solutions  are  found  to  the  linearized  three  dimensional  equations 
for  free  surface  gravity  waves  trapped  against  a  straight  coastline  with  a 
variable  (perpendicular  to  the  coastline)  topography.  Three  families  of 
topographies  are  found,  one  concave  upwards  and  two  convex  upwards,  which  will 
support  these  edge  waves  as  separable  solutions  to  the  original  equations. 

For  a  given  topography,  specified  by  an  initial  slope,  h^»  a  typical  water 
depth,  h^,  and  a  typical  offshore  distance,  x^,  solutions/ are  given  in  terms  of 
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the  nondimensiooal  parameters,  h^  x^/h^',  hV  *The  first  parameter  is  a 


/Chat 


measure  of  how  the  initial  slope  h*  compares  tO/tnat  for  a  straight 
topography  from  the  origin  to  (x^,h^),  i.e.  ^(h^/x^) .  This  parameter  characterizes 
t.he  families  in  ranges  0  <  h^  <  1/2,  1/2  <  h^  x^/h^  <  1  ant* 

1  <  h^  Xjj/h^..<  (in  10)/0.9.  (This  latter  constant  bound  is  a  function  of  the 
definition  of  offshore  scale  and  can  be  modified  to  other  values.) 

'™?The  nondimensional  Crequcncv,  period,  vertical  wave  number,  offshore  decay  rate 

'  jr*  '  $  u.  If  ft 

and  topography  can  be  expressed  in  terms  cf  the  single  parameter  h^  x^/h^; 

jrxi  o  ** 

the  velocities,  wave  height  and  alongshore  wave  number  depend  on  both  parameters. 
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INTRODUCTION 

The  trapping  of  wave  energy  along  coastlines  has  been  of  theoretical 

V 

interest  since  the  original  edge  wave  solution  of  Stokes  (1846). 

Recently,  there  has  been  consideration  of  the  physical  reality  of  such  trapped 
waves  and  a  number  of  observations  have  led  to  a  belief  that  such  trapped  waves 
may  indeed  play  a  role  in  beach  erosion,  local  resonances,  anisostatic  response 
to  traveling  pressure  disturbances,  etc.;  a  review  of  these  pcints  is  contained 
in  LeBlond  and  Mysak  (1977)  and  (1978).  Some  distinction  must  be  made  at  the 
outset  of  any  research  as  to  the  type  of  trapped  waves  to  be  discussed. 

Here  the  emphasis  is  on  class  1  waves,  i.e.  gravity  dominated,  of  frequencies 
high  enough  that  Coriolis  effects  may  be  neglected,  i.e.  edge  waves. 
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These  waves  arc  basically  high  frequency,  dispersive  and  can  travel  in  either 
direction  along  a  coastline;  in  fact,  since  no  rotational  effects  are  included 
and  thus  no  direction  is  preferred,  the  coast  line  may  be  oriented  in  any 
direction. 

The  problem  at  hand  then  is  to  examine  the  existence  of  edge  waves 
traveling  along  and  trapped  against  a  straight  coastline  bounding  a  semi-infinite 
ocean  whose  depth  contours  arc  parallel  to  the  coastline.  The  ocean  is 
considered  to  be  a  perfect  homogeneous  fluid,  inviscid  and  irrotational. 

The  purpose  of  the  paper  is  to  establish  exact  solutions  for  these  three 
dimensional  equations  by  reversing  the  usual  question  "what  are  the  trapped 
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solutions  for  a  given  topography"  to  read  "what  topographies  will  support  trapped  J 
wave  solutions".  One  solution  already  available  for  comparison  is  that  of  :1 
Stokes  (1846)  who  considered  a  linear  topography  (constant  slope)  aad  obtained  a  1 

'  1  -ffl 

"zeroth"  mode  solution  which  Ursell  (1952)  extended  to  a  number  of  discrete  nodes  I 
plus  a  continuous  spectrum.  This  was  based  on  the  previous  work  of  Eckart  (1951)  J 
who  obtained  higher  modes  using  a  shallow  water  (2D)  approximation. 


Other  authors  [e.g.  Reid  (1958)]  added  Coriolis  effects  to  the  shallow  water 
theory.  In  fact,  the  bulk  of  the  effort  in  this  area  has  relied  on  the  shallow 
water  assumption  [e.g.  see  p.  460  LeBlond  and  Mysak  (1977)]  and  the  lack  of 
other  exact  three  dimensional  edge  wave  solutions  was  commented  upon  by 

Grimshaw  (1974)  who  examined  upper  and  lower  bounds  on  such  dispersion  relation- 

0: 

ships.  Some  non-linear  topographies  have  been  considered.  Robinson  (1964) 
examined  class  11,  i.e.  quasi-geostrophic  low  frequency,  waves  on  a  linear 
sloping  shelf  of  finite  width  terminated  by  a  discontinuous  drop  to  a  constant 
depth  infinite  ocean.  Ball  (1966)  used  an  exponential  depth  profile  with 
solutions  in  the  form  of  hypergeometric  or  Jacobi  polynomials  to  examine  both 
edge  and  continental  shelf  waves,  i.e.  both  class  1  and  class  II  waves. 

Mysak  (1968)  also  discussed  both  classes,  using  Robinson's  depth  profile,  where 
the  solution  is  given  in  terms  of  Laguerre  polynomials.  Hidaka  (1976-a) 
examined  seiching  due  to  a  submarine  bank  described  by  a  parabolic  depth 
variation  in  terms  of  Mathleu  functions.  This  is  a  trapped  but  not  specifically 
an  "edge"  wave,  but  it  is  of  class  I  nevertheless.  In  a  related  paper 
(Hidaka  (1976-b))  he  examined  shelf  resonances  when  the  apex  of  the  parabola 
reaches  the  free  surface.  Murty  (1977)  gives  a  general  review  of  tsunami 
theory  and  current  research,  which  includes  theoretical  and  observational 
studies  on  edge  waves  caused  by  tsunamis  (which  are  generally  of  high  enough 
frequency  that  Coriolis  effects  may  be  neglected,  i.e.  class  I  waves)  but  are 
long  enough  that  a  two  dimensional  shallow  water  theory  suffices. 

Experiments  performed  by  Ursell  (1952)  indicated  the  reality  of  such  edge 
waves;  further  experiments  by  Galvin  (1965)  indicated  that  these  could  be. 
excited  (nonlinearly)  by  normally  incident  waves.  Huntley  and  Bowen  (1973)  and 
Huntley  (1976)  have  made  observations  which  indicate  the  presence  of  edge  waves 
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along  coastlines  in  the  U.K.  while  Nakamura  (1962),  Hatori  and  Takahasi  (1964) , 
Nakamura  and  Watanabe  (1966),  Aida  (1967),  Hatori  (1967)  and  others  have 

t- 

described  the  excitation  of  such  edge  waves  by  tsunamis  incident  onto  the  coast 
of  Japan.  There  is  other  observational  evidence,  Wilson  and  Toruo  (196R),  that 
some  tsunami  energy  is  trapped  on  the  continental  shelf  in  the  generating  region 
at  least  for  the  1964  Alaskan  Earthquake  and  it  may  be  expected  that  some  of 
this  energy  could  be  converted  to  and  travel  along  the  coastline  as  edge  waves. 
Olsen  and  Hwang  (1974)  indicate  that  edge  waves  may  play  a  significant  role  in 
near  shore  tsunami  behavior. 
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FORMULATION 


The  governing  equations  are  veil  known,  e.g.  Lamb  (1932)  and  are  written 
here  for  a  straight  coastline,  x  ■  0,  -«  <  y  <  4®,  bounding  a  semirinfinite 
ocean,  x  >  0  with  a  free  surface  at  z  ■  0  and  a  bottom  topography,  z  «  -h(x) 
with  *  positive  upward: 

(1)  V2$(x,y,z,t)  -  0;  x  >  0,  -h  <  *  <  0,  -®  <  y  <  » 

(2)  o2$  -  g  3*/3z  -0;  x  >  0,  z  ■  0,  -•  <  y  <  • 

(3)  3$/3 z  *  -(dh/dx)  3$/3x;  x  >  0,  z  *  -h,  -•  <  y  <  «* 

where  a  time  harmonic  behavior,  exp(-iot)  has  been  assumed  for  the  velocity 
potential,  The  Coriolis  effect  has  been  neglected  and  the  y  dependence  will 
be  taken  as  periodic  (the  edge  wave  assumption) .  Assuming  separation  of 
variables,  the  velocity  potential  is  given  as  $  ■  X(x)  Z(z)  cos  (ky-ot) ,  using 
a  real  form  for  the  and  C  dependence,  leading  to: 

(4)  X(x)  -  A  exp  (-lx)  +  B  exp  (+ix) 

(5a)  Z(z)  »  C  sinh  (nz)  +  D  cosh  (nz) 

(5b)  •  C*  sin  (mz)  +  D’  cos  (mz) 

where  t2  is  a  separation  constant  and  n2  -  k2  -  l2  >  0  and  m2  ■  i2  -  k2  >  0 

provide  two  different  solution  forms  for  Z(z).  The  coefficient  B  is  taken  to 
be  zero  to  provide  a  "trapped"  solution  and  A  may  be  absorbed  into  Z(z). 

Substituting  the  first  solution,  (5a),  into  the  boundary  conditions  provides 

(6)  C/D  -  o2/gn 

(7)  dh/dx  -  -(n/i)  -teak aVirL. 

Kn  h/  WK}  [1  -  (oi/gn)tanh(nh)] 


(«) 

(b) 


-(n/t)tanh(nh-tanh  ^(o2/gn) ] ;  o2/gn  <  1 
-(n/l)coth[nh  -  coth  ^(c2/gn)];  o2/gn  >  1 
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Using  the  boundary  conditions  in  (5b)  leads  to 

(8)  C'/D'  »  o2/gm 

(9)  dh/dx  «  (n/i)tanrmh  +■  tan  ^(o2/gm)] 

Equations  (7)  and  (9)  can  be  Integrated  to  give  three  families  of  topographies 
that  give  trapped  solutions  as  summarised  below: 

Case  I:  k2  >  i2;  o2/gn  <  1 

Case  II:  k2  >  l2;  o2/gn  >  1 

Case  III:  k2  <  i2;  all  o2/gm 
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SPECIFIC  SOLUTIONS 

Case  I  leads  to  a  concave  upward  topography  with  h(0)  set  equal  to  zero 
(to  prevent  flux  across  the  coastline) ,  an  asymptotic  constant  depth, 
hw  ■  (tanh  *(o2/gn)]/n  and  an  initial  slope,  h^  s  (dh/dx)Q  »  a2/gt. 

Using  0  ■  tanh  *(o2/gn) , 

h(x)  ■  [0  -  sinh  *(sinh  0  •  exp(-n2x/f.])  )/n 
with  a  velocity  potential  given  as 

$  »  D  exp  (-lx)  cos  (ky-at)cosh[n(z+ha|)  J/coshCnh^) 

Case  II  (k2  >  £2,  o2/gn  >  1,  i.e.  eq.  7b)  leads  to  a  convex  upward  topography  with 
an  infinite  slope  at  a  finite  value  of  h  and  x.  Using  0'  *  coth  *(o2/gn), 
h(x)  *  (O'-cosh  ^(cosh  0* *exp[-n2x/i])]/n 
The  maximum  extent  of  the  topography  is  Xj^  -  (t/n2)in[cosh  O’]  which  occurs  at 
h^  ■  0'/n.  The  velocitv  potential  in  this  case  is,  for  x  <  x^, 

$  ■  D  exp(-ix)cos(ky-ot)sinh(n(z+hM)]/sinh(nhM) 

When  o2/gn  «  1,  both  cases  I  and  II  recover  the  Stokes  solution,  with 
h  ■  x  tan  o,  l  ■  -k  cos  a,  n  *  -k  sin  a  and  a  »  tan  ^(h^). 

Case  III  (k2  <  l2,  all  o2/gm;  i.e.  eq.  9)  using  y  -  tan  1  —  leads  again 

o® 

again  to  a  convex  upward  topography  with  an  infinite  slope  at  «  -(i/ra2)fn[siny] 
with  depth  h^  *  (*/2  -  y)/m 

h(x)  *  (sin*’i(sin  y  •  exp(4w2x/£])-y]/m 
with  a  velocity  potential  for  x  <  given  by 

$  ■  D'  exp(-ix)cos(ky-ot)sin[m(h}il+x))/sin(mhM) 

These  last  two  families  of  solutions  both  appear  to  be  physically 
unrealistic  in  that  the  prescribed  bottom  topography  does  not  extend  past  a  finite 
distance  and  finite  depth  before  "doubling  back."  If,  however,  the  resulting 
velocity  fields  for  these  problems  ljave  decayed  sufficiently  with  distance  from 
the  straight  coastline  through  the  term  exp(-ix)  and/or  if  they  have  decayed 


£§§3$i^lW!l58!f£*^^  M^-c-  ~, 
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sufficiently  with  depth  to  have  both  u  and  w  negligible  at  the  lower  boundary 
before  these  "turning  points"  are  reached,  the  point  is  moot  since  the  lower 
boundary  condition,  cq.  (3),  is  essentially  satisfied  identically  for  all  x  and 
z  past  these  turning  points.  Thus  these  solutions  may  provide  a  basis  for 
realistic  trapped  waves  on  a  topography  which  matches  that  prescribed  by  these 
solutions  prior  to  the  turning  points  and  a  flat  (or  any  other)  topography  past 
them. 

The  rest  of  this  report  will  examine  the  implications  of  these  solutions, 
particularly  for  the  various  limiting  cases.  To  do  so  conveniently,  it  is  useful 
to  nondimeiisionalize  the  parameters  and  variables  used.  The  coordinates  x,  y 
and  z  are  nondimtnsionalized  with  respect  to  h^Ch^  **  h^  for  case  I)  as  are  the 
parameters  k,  i.  and  n  (a  will  be  counted  with  n  generically) , 
i.e.  (x,y,Z)  ■  (x.y.zj/h^  and  (k,i.,n)  -  (k,£,n)  (h^)  as  the  case  requires. 


The  frequency  is  scaled  to  oq  ■  (g/h) 


1/2 


as  is  the  time,  i.e.  5  *  a/aQ  and 


C*  ■  t  ♦  o  .  Velocities,  V  *  V$,  are  scaled  by  u  *  u(x"0,  y*0,  z*0,  t"0)  and 


the  wave  height,  £  -  -(l/g)|^|z  ■  0),by  C  -  C(x»0,  y-w/2k,  z«0,  t-0)  which 
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equals  (o/gt)u  • 
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SOLUTIONS  AND  LIMITING  CASES 

The  above  solutions  ere  defined  by  three  dimensional  parameters  or  two 
nondiaensional  parameters,  using  h^  as  a  typical  length  scale.  The  choice  of 
those  parameters  to  be  prescribed  is  important  since  it  establishes  how  simply  the 
results  may  be  described.  One  such  nondiaensional  "fundamental"  parameter  will 
be  chosen  as  the  initial  slope,  h^.  The  second  choice  depends  on  which  aspect 
of  the  problem  is  under  study.  For  examination  of  various  "asymptotic  solutions," 
the  parameter  3 2/fl  (or  0,  9'  or  y  )  arises  naturally  in  the  derivation  and  will 
therefore  be  used  for  these  cases.  The  remaining  computations  for  the  general 
solutions  however  will  be  carried  out  using  h^  x^  as  the  second  "fundamental" 
parameter.  For  case  I,  is  defined  as  h(x  ■  .9  x^)  *  .?  hw;  actually  other 
choices  could  equally  well  be  made,  to  define  a  typical  offshore  length,  e.g. 
h(x  -  (1-5  )x^)  •  (1-5  )hw.  Then  is  x^/h^  and  the  parameter  h^  is  the 

ratio  of  the  actual  Initial  slope  to  the  slope  of  a  straight  line  from  the  origin 
to  (x^,h^).  This  parameter  characterizes  the  three  families  as  well  as  did 
32/n. 

Case  I  requires  1  <  h^  x^  <  (in  10)/.9  (or  more  generally,  (in  l/5)/(l-5)) 
which  is  a  curve  whose  initial  slope  h^  is  greater  than  hj,/^  giving  a  curve 
lying  below  the  straight  line  from  the  origin  to  (.9  x^, 
i.e.  is  concave  upwards.  Case  II  and  case  III  require  0.5  <  h^  ^  <  1  and 
0  <  h^  x^  <  0.5  respectively;  these  both  give  topographies  lying  above  (closer 
to  the  surface)  that  the  straight  line  to  are  concave  downwards. 

Thus  the  three  families  can  be  represented  by  a  continuous  range  of  values 
of  h’  x^,  with  some  question  as  to  the  cross  over  values  of  1.0  (which 
corresponds  to  the  Stokes  solution)  and  0.5  (which  presents  no  actual 
difficulty).  The  main  advantage  to  use  of  this  parameter  is  that  the  frequency. 


.9  h^)  for  X  <  x^. 
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U,  the  vertical  "wavenumber,"  5,  the  (modified)  offshore  decay  rate,  T/x^,  *' 
and  the  topography,  HCx/ic^)  can  be  represented  as  functions  of  h’  ^  alone; 
i.e.  they  are  one  parameter  functions.  The  along' shore  wavenumber,  k,  the 

t  .  f  aV 

wave  height,  f,  and  the  velocities,  u,  v,  w,  depend  oh  both  h^  and 

h^,  i.e.  are  two  parameter  functions.  Thus  much  of  the  solution  cab  be  expressed 

quite  simply  in  terms  of  a  single  parameter. 


-./■a 
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CASE  X: 

This  is  in  many  ways  the  most  satisfying  solution  with  «  well  behaved 
topography.  Consider  h^,  h^  and  x^  ■  (l/.9)x(h  »  *91^)  to  be  specified. 

In  non-dimensional  terns,  R  is  given  by  and  since  is  1  by  definition: 


tanh  n 


in 


fsinh  5 


sinh‘9R 


'  •’  f* 


with  3,  t  and  k  given  by 


3  ■  [n  tanh  n]^2 
l  -  [n  tanh  n]/h^ 
fc  -  R(1  +  (tanh2fi)/h^2]1/2 


If  o2/n  were  specified  instead  of  X^,  R  would  be  s.laply 
tanh“1(52/n)>  eliminating  the  need  to  solve  a  transcendental  equation  for  n. 
The  topography  is  given  by 

h{x)  -  1  -  (l/n)sinh_1(sinh  R  •  expC-n2*/!)) 


and  the  velocities  and  wave  height  oy 

u  *  exp(-tX)cos(k?  -  3t)coshtR(l+i)]/cosh  R 
?  -  (k/i)exp(-Ix)sin(ky  -  3t)cosh[n(l  +  E)]/cosh  n 
w  -  -(n/t)exp(-ix)cos(ky  -  3t)»inh[R(l  +  z)]/cosh  n 
X  "  exp  (-lx)  sin  (icy  -  ot) 

( 

Note  chat,  as  5  approaches  zero,  h^  5^  approaches  1  while  as  n  approaches 
infinity,  fcj  5^  s^pregshes  (In  10)/. 9.  Furthermore,  R  and  3  depend  only  on 
the  combination  h^  5^,  while  I,  k,  h  and  the  velocities  require  separate 
values  for  h^  and  x^.  However,  l/x^  and  hix/x^)  depend  only  on  h^  x^* 


A:  Consider  the  limiting  case  32/n  *  e  «  1.  Approximation*  for  the  governing 
parameters  are  found  to  be 

5«e 

I  .i*  c2/h' 
o 

5  C 

k  *c(l  +  (c/h')2]l/2 
o 

h'  Xy,  *  (1/.9)  {in  10  **  c2((ln  10)/3  -  1/  )] 
on 

The  restriction  on  o2/n  implies  o  <<  aQ  ,  i.e.  that  a  shallow  water  wave 
theory  be  used  leading  to  the  same  topography  for  all  applicable  frequencies. 

The  topography  is  given  by 

h(x)  -  1  -  exp[-V  *i 

which  in  turn  is  the  same  ns  that  used  by  Ball  (1967).  Within  the  limitations  of 
this  case,  the  velocity  potential  is 

$  •  D  exp(-K?X/h^)cos(*:(l  +  (>c/h^2)^2y  -  ot) 

with  k  B  8 .  Tills  corresponds  to  Ball's  solution  for  the  case  of  sero  rotation 
and  node  number  (n)  *  0,  i.e.  the  first  term  in  hie  hypergeometric  polynomial 
which  then  reduces  to  a  simple  exponential  decay  in  x.  The  corresponding 
velocity  field  has  no  vertical  variation,  an  exponential  decay  with  distance 
perpendicular  to  the  shore  line  and  a  periodic  wave  behavior  along  the  shore 
line,  both  of  which  depend  on  the  nondimens lbnal  frequency  and1  the  initial 
beach  slope.  Since  this  form  has  already  been  obtained  by  ball,  in  far  greater 
generality,  it  will  cot  be  discussed  further  here.  ~  ^ 
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B:  Consider  next  the  case  when  32/n  is  near  1,  i.e.  a2/ n  »  1-e;  c  «  1. 

Then  H  *  ^tn(2/e)  Is  a  large  nuaber.  The  remaining  parameters  follow  as 

o  =  [|£n(2/c)]1/2 

I  ^  (1/2  h^)tn(2/c) 
k  •*>  (1/2)  £n(2/ c)  [1  +  (l/h^)1/2 

and 

h^  JCjj  o»  1  -  ln(l  -  (c/2) *1)/(.45  in(2/c)) 
i.c.  h^  is  slightly  larger  than  lv 

The  topography  is  given  by 

(h’R-l)/2“| 

F.(x)  1  -  (l/2tn(2/c))  ^ainh  -L[jCc/2)  °  J 

For  h^X  «  1,  this  reduces  to  a  linear  topography,  fi  *  h^x. 

Since  h’0  is  close  to  l,  this  implies  that  x  «  x^  for  this  linear 

topography.  For  these  small  values  of  x,  the  velocities  and  wave  height  are 

3  *  (2/c)1/2(i"X/ho)cos(y«(l/2tn(2/c)(l+(l/h*)2)1/2  -  (l/2tn(2/e))1/2  t) 

o 

5  *  (1  +  h^2)(2/c)1/2(*~S/h^  sin  1  "  "  J 

w  ^.-h,(2/c)1/2(Z_5/hi)  cos  t  "  "  J 

o 

l  -  (2/c)"x/2h°  sin  f  ”  "  ) 

This  solution  recovers  the  classical  Stokes'  solution,  but  only  for  large 
n  and  3.  However,  if  the  general  solution  is  examined,  the  complete  Stokes' 
solution  (for  all  3)  is  obtained  for  S2/B  exactly  equal  to  1.0.  This  limiting 
process  for  case  I  can  only  provide  the  high  frequency  Stokes' solution  since  the 
regular  Stokes'  solution  is  actually  neither  in  class  I  nor  in  class  II  but  on 
their  boundary.  Since  1  <  0,  these  velocities  di/a  out  rapidly  with  increasing 
depth  and  distance  offshore. 
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CASE  II: 

This,  together  with  case  III,  compose  the  solutions  for  h’  <  1, 
i.e.  the  concave  downward  topographies;  n  is  related  to  through 

h;  %  -  ?■  .  in(cosh  n) 

* 

with 

-1/2  * 

5  -  (ft  coth  n) 

l  •  (B  coth  n)/h^ 

k  ■  5(1  +  (coth2S)/h'2]*^ 

o 

and 

h(5)  -  1  -  (l/n)cosh-*[cosh  n  •  exp(-tl22/Z) ] 

The  velocities  are  then  given  by 

u  -  exp  (-lx)  cos  (Icy  -  5t)sinh(n(l  +  z)]/sinh  n 
?  ■  (k/t)exp(-lx)ain(5cy  -  ot)sinh[n(l  +  z))/ainh  n 
w  -  -(n/Z)exp(-Ix)cos(Icy  -  ot)cosh[n(l  +  z)]/sinh  n 

and  the  wave  height  by 

(  *  exp  (-lx)  sin  (icy  -  at) 

As  n  approaches  zero,  hSc^  approaches  1/2  while  as  n  approaches 

infinity,  h^  approaches  1  from  below.  Again, 

depend  only  on  the  single  parameter  h'  xu. 

o  n 

Limiting  values  of  o2/n  may  again  be  examined. 

A:  Consider  the  limiting  case  32/n  •  1  +  c  where  c  «  1.  Here,  n  l/2jtn(2/e) 
ag.ln  and  ^  ^  1  '  f^i7cT  vlth 


5,  i/Xjj,  hCx/j^)  and 


and, 
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o  «  (l/21nC2/c) l1^2 

I  *  (£n(2/E))/2h; 

1/2 

k  s»  (l/2h^)£n(2/e) (1  +  h;2)A/ 

for  xh'  <  1  (aa  it  must  be  for  x  <  3c.) 
o  n 


h(x)  **•  x 

u  of  (2/c) ~  K/ho)/2  cos(y*(l/2  in(2/c)*(l  + 

d/h;)2)1/2) 

v  -  (1  +  h;2)1/2(2/c)(s  ‘  >:/hi)/2  aint 

n 

V  :  -h^(2/C)(Z  "  R/h«)/2  COS  ( 

n 

l  ,  (2/c)~x/2h°  sin  ( 

m 

(1/2  *n(2/e))l/2  t] 
"  1 

"  ] 

M  1 


For  small  x,  3,  the  Stokes  solution  for  large  3  ia  found  again  as  would  be 
expected.  In  fact,  this  solution  is  identical  to  1-B  except  for  ^  which  is 
greater  than  1  for  I-B  and  less  than  1  here  for  non-xero  c. 
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B:  Next  consider  the  limiting  case  o2/n  •  1/c  where  e  «  1.  Here,  n  is 
approximately  c  +  (1/3) c3,  h’  x^  1/2  +  e2/12  and 

5*1  +  e2/6 
I  (1/h^Xl  +  c2/3) 
k  [c2  +  (l/h^)2]1/2(l+c2/3) 

and 

h(X)  1  -  (1  -  2h*  x)1/2 
o 

which  for  small  h©  JC  is  linear  again  with  slope  hX  The  velocity  field  and 
wave  height  are  given  by 

a  (1  +  E)cxp(-x/h^)cos(|7*(e2  +  (i/h^)2)1^2  -  t] 

v  ~(1  +  (e  h')2]1/2(l  +  2)exp(-Jt/h’)ain[  «  ] 

o  o 

w  -  h' •exp(-x/h,)cos{  "  ] 

o  o 

l  exp(-x/h^)sin(  "  ] 

Note  that,  to  the  lowest  order  in  c,  h^  R^,  o  and  I  are  essentially 

independent  of  c,  as  arc  k,  u,  v  and  w  if  h^  «  1  as  well. 

The  topography  is  parabolic,  with  a  maximum  horizontal  extent  of  1/2  as 
expected.  The  horizontal  velocities  decrease  linearly  with  depth,  but  do  not 
vanish  at  the  bottom  until  h  equals  h^,  i.e.  z  can  reach  -1. 

The  vertical  velocity  is  constant  with  depth  and  is  0(hg)  compared  to  the 
horizontal  velocities;  this  is  then  a  small  value  unless  x^  is  itself  of 
order  0(1),  i.e,  the  initial  slope  is  45*  or  greater  which  is  physically 


unlikely. 


CASE  III:  (k2  <  i2) 


Here  the  range  of  x^  it  from  0,  as  m  approaches  ir/2  to  1/2,  as 
5  approaches  zero  with  the  general  relation 


cot  5«in(cos  5) 


-  h*  5L 


o  -  (5  cot  ra)^2 


£  *  (m  cot  in)/h' 
c 


s[m2  - 1 


sin  ^fsln  yexc 


Y  •  -5  +  i?/2  and 


u  -  expC-Ix) *coB(ky  -  ot) *sin(5(l  +  z)]/sin  ra 


v  -  (K/I)cxp(-Ix)sin(ky  -  5c)sin{m(l  +  z)]/sin  ra 


w  ■  -  (5/ 1)  exp  (HU)  cos  (ky  -  ot)cos [5(1  +  z)]/sin  ra 


£  -  exp(-Ix)sin(ky  -  ot) 


A:  The  first  limiting  case  here  is  o2/I  *  c  «  1.  Then  y  ~  g  and 


5  x  n/2  -  c  with  ~  - 


2c  in  c 


o  W  (n  e/2)J 


I  *  t  b/2  h* 
o 


k  *(*/2)[(e/h')2  -  1]1/2 


gaafWhi. 
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ii(x)  u  2/r(sin  i(exp[-  n(x^  -  x)h^/2ej)  -  c] 


Clearly  e  >  h*  in  order  to  have  real  k,  i.e.  a  progressive  wave  solution, 
o 


If  x  is  snail  enough  such  that  exp(axh^/2e)  1  +  *8h^/2c  the  linear  topography 


solution  is  round  again,  i.e. 


h  ar  h'  X 
o 


The  velocities  and  wave  height  are 


u  a  cxp(-etrx/2h!) *cos[y*r/2*((e/h')2  -  l)^2  -  (ire/2)^2  t]sin[7r/2(l  +  £)] 


v  i  [1  -  (h,/c)2]1/2*exp(-effif/2h,)sin[ 
0  0 


]sinfir/2(l  +  *)] 


w  «r  -(h^/c)exp(-c5Jx/2h^)cos( 


]cos[w/2(l  +  I)] 


1  ~  exp(-crrx/2h^)sin( 


Such  a  solution  requires  extremely  snail  values  of  the  initial  slope,  h'. 


This  in  turn  leads  to  a  rapid  decoy  of  wave  heights  and  velocities  in  the  offshore 


direction  due  to  the  exp[-crx/2h^] ,  e.g.  for  e  equal  to  h^  (such  that  k  is 


zero  corresponding  to  no  alongshore  variation) ,  the  velocities  have  decayed  to  AX 


of  thier  shoreline  values  when  x  has  reached  2  h^.  The  wave  periods  are  very 


long,  which  could  cause  Coriolis  effects  to  become  significant  and  thus 


invalidate  this  solution. 


B:  The  final  limiting  case,  c2/fi  *  1/e,  e  <<  1  implies  y  a.  if/ 2 -e,  fi  ■fc  e  and 
».»  s  ~  i /•>  ..ui.  r:  ..  i  ~  t i/u*n  C  ~ 


h'  5L,  ts  1/2  with  o  *  1,  I  *  (l/hf),  It  c  (l-(eh’)2)  '  /h*  and 
on  o  o  o 

h(x)  ~  1  -  (1  -  .‘^h’}1^2  which  reduces  to  the  linear  form  when  x  h’  is  email. 


The  velocities  and  wave  height  in  this  case  are 


mm’W^sszm;?* 
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Q  «  a  +  *)«xp(-*/h’)cos[y<a  -  (eh»)2)1/2/h.’  -  t] 

O  0  0 

v  Ml  -  (ch^)2)“^2(l  +  z)exp(-x/h^)sin[  "  ] 

w  ss  -h^*exp(-5c/h^)cos[  "  I 

;,  «*  exp(-x/h*)sin[  "  ) 

which  match  the  limiting  case  II-B  as  they  must  since  both  of  these  correspond 
to  the  case  h^  ■  1/2  which  is  common  to  both  case  II  and  case  III. 

Apart  from  the  alongshore  velocity,  7,  and  the  corresponding  alongshore  dependence 
on  ky,  this  solution  is  independent  of  e.  Furthermore,  as  c  passes  through 
zero,  the  match  between  case  II-B  and  case  III-B  is  complete,  i.e.  the  transition 
from  case  II  to  case  III  is  smooth. 
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DISCUSSION  OF  GENERAL  SOLUTION 

Although  Che  problem  is  completely  specified  by  three  dimensional 
(e.g  h^,  ^  and  h^)  or  two  nondimensional  (e.g.  h^  and  x^),  it  is  of 
practical  interest  to  express  results  in  terms  of  as  few  parameters  as 
possible.  This  may  be  done  for  a  portion  of  the  present  problem. 

The  quantities  R,  5,  I/x^  are  functions  only  of  the  single  parameter 
(h^  x^)  and  the  topography  ii  is  a  function  of  the  parameter  h^  x^  and 
the  spatial  variable  x/x^.  Thus  although  the  computer  program,  given  as 
appendix  A,  uses  three  dimensional  inputs,  h^,  h^  and  x^,  the  results  may 
be  stated  in  a  fairly  compact  form. 

The  vertical  "wavenumber,"  R,  is  related  to  the  single  parameter 
h^  *M  ^  h'0  “  F(n)  where  F(R)  has  different  forms  for  the  three  cases 

considered  and  is  plotted  in  Figure  (1).  Thus  n  may  be  read  from  this 
figure  for  any  value  of  h^  between  0  and  (in  10)/0.9 
(or  calculated  through  the  program).  Figure  (2)  provides  the  (nondimensional) 
period,  T,  as  a  function  of  the  single  parameter,  h^  Jt^.  Since  the 
nondimensionalization  involves  h^,  dimensional  periods  require  knowledge  of 
both  h^  and  h^  x^.  At  h^  x^  equal  to  1,  there  is  an  anomaly  in  the 
original  solution.  At  this  value,  equation  (7)  reduces  to  dh/dx  equal  to 
a  constant  (-n/l)  and  h^  is  1  for  all  values  of  h^.  In  a  sense,  this 
ia  a  "fourth"  solution  family,  the  well  known  Stokes  solution,  in  that  it 
can  not  be  obtained  from  the  other  solutions,  except  as  a  high  frequency 
limit.  Thus  Figure  (2)  has  zero  period  at  h’  equal  to  1  with  the 
Stokes  solution  indicated  as  a  dashed  line.  The  Stokes  solution, 
h  •  +  tan  a  •  x,  requires  o2  *  gk  sin  a.  Thus  for  any  h^  *  tan  o,  there 
is  a  single  specific  frequency,  o,  at  least  for  the  lowest  separable  mode. 
Ursell's  extension  (1952)  to  higher  modes  is  not  comparable  here  since 
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those  higher  modes  were  not  separable  solutions  and  could  not  be  obtained 
by  this  present  approach. 

A  question  then  arises  as  to  the  appropriate  solution  for  an  "almost" 
linear  topography.  The  examples  here,  which  represent  the  only  separable 
solutions— but  clearly  not  necessarily  the  only  solutions — will  only 
support  very  high  frequency  trapped  waves  when  h^  2^  is  near  to  1  but  not 
equal  to  1. 

The  offshore  decay  parameter,  i,  depends  on  both  h^  x^  and. 
but  in  a  manner  such  that  i/x^  is  a  function  of  h^  2^  only.  This  is 
shown  in  Figure  (3)  where  I/x^  is  given  for  h^  -  200  meters  atwl 
Xjj  *  200  kilometers  for  various  values  of  h^  but  valid  of  course  for  any 
choices  of  and  x^.  This  figure  indicates  the  limited  range  of 
practical  solutions  since  all  terms  decay  like  exp(-Ifc)  or  equivalently 
expt-d/Xjj)  (x/x^Xjj2) .  To  have  this  exponential  greater  than  0.04 
compared  to  a  value  of  1  at  x  *  0,  (I/x^)  <  ir/Cx/x^Hx^)2.  Even  for 
x/x^  up  to  0.1  representing  the  region  of  measurable  wave  heights,. 

(I/Xjj)  <  31.4/Xjj2.  Realistic  choices  for  5^  are  of  the  order  of  .40  to 
1000  for  the  continental  slope  and  shelf  respectively;  these  limit 
£/xjj  <  2  x  10  and  3  x  10  which  in  turn  imply  that  h^  X^.must  be  close 
to  either  zero  or  (in  10) /0. 9  in.  either  case.  Thus  for  realistic  ocean 
topographic  scales,  the  energy  oust  be  trapped  very  near  shore. 

Figure  (4)  shows  the  topography,  h,  as  a  function  of  2/x^  for 
several  values  of  the  parameter  h^  x^»  based  on  a  choice  of  h^  equal 
to  200  meters  and  equal  to  200  kilometers  but  again, valid,  in,  general. 
Those  topographies  for  >  i  are  forced  to  pass,  through  (0.9,  0r9)  by 

the  choice  of  definition  for  x^;  other  choices  would  simply:  rescale  these 
curves  to  pass  through  (1-5,  1-5),  i.e.  appear  steeper  for  5  <  0.1. or, 
shallower  for  5  >  0.1.  As  V  x^  approaches  1.0,.  the  topography  appears 
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almost  linear  except  near  x/x^  equal  to  1.  For  values  of  55^  less 

than  1,  convex  upward  topographies  are  found  with  an  infinite  vertical 

slope  at  x/5^  equal  to  1.0.  The  transition  across  h^  equal  to  0.5 

is  smooth  with  no  ascertainable  peculiarities. 

Calculations  of  velocities  and  waveheights  produced  extremely  small 

values  for  any  offshore  distance  that  itself  was  not  a  very  small  fraction 

of  For  example,  with  h^  ■  0.002556,  h^  *  200  meters,  * 

x^  •  200  kilometers  the  velocity  and  waveheight  values  had  reached  4%  of 

their  maximum  values  when  x  reached  0.2%  of  x^,,  e.g.  about  400  meters— 

at  which  distance  the  water  depth  was  approximately  one  meter.  This  implies 

that  the  edge  wave  was  completely  trapped  in  the  region  where  the  model  is 

the  poorest  representation  of  the  overall  topography  and  is  therefore 

Inadequate.  For  this  reason,  no  results  are  shown  for  this  case. 

For  much  smaller  values  of  however,  it  is  possible  to 

illustrate  the  depth  dependence  for  velocities.  Figures  (5)  and  (6)f show 

results  for  x^  -  1.0.  The  first  plots  the  offshore  velocity,  u,  as  a 

function  of  (.'./x^)  and  z  for  h^  *  1.01.  The  second  illustrates  the 

different  coastline  values  for  v  and  w  as  compared  to  the  reference 

value  of  1  for  u,  as  a  function  of  h'  5L..  There  is  a  cutoff  at 

on  y 

h^  -  0.2854  for  x^  ■  1.0  below  which  k  becomes  imaginary.  This  corre- 

sponds  to  the  requirement  that  m  >  !n(see  15).  Similarly  near 

h'  iL,  -  (la  10)/0.9,  V  behaves  like  (1  +  (h’/c)2)1/Z,  as  defined  in  case  I-A 
on  o  *  ,  • 

while  w  behaves  like  -h^. 
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CONCLUSION 

The  bulk  of  observations  have  been  for  shallow  water-long  wavelength 
situations  such  Chat  the  question  of  vertical  variation  does  not  enter,  as 
befits  the  general  emphasis  on  two  dimensional  theory  for  edge  waves. 

All  tsunami  observations  are  of  this  form  as  were  the  observations  of  Huntley 
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and  Bowen  (1973)  who  found  evidence  for  the  lowest  mode  edge  wave  of  Ball’s 
(1967)  solution  corresponding  to  case  I-A  here.  Ursell  (1952)  described 
experiments  on  a  linear  topography  but  did  not  give  any  details  of  vertical 
structure.  Thus  it  appears  that  separate  observations  and/or  experiments 
may  be  necessary  to  detect  the  three  dimensionality  of  those  edge  waves 
described  here.  Whether  such  observations  are  warranted  is  another  question. 
Wavelengths  that  are  significant  enough  to  carry  appreciable  energy 
alongshore  may  also  need  to  be  long  enough  to  allow  a  two  dimensional  theory 
to  apply.  Such  a  theory  allows  for  far  more  general  topographies  than  does 
the  three  dimensional  theory,  as  discussed  in  appendix  B.  This  is  clearly 
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a  question  for  further  consideration  on  physical  grounds. 

Assuming  this  research  to  be  physically  justified,  there  are  still 
several  mathematical  questions  left  unanswered.  First,  are  there 
non-separable  solutions,  analogous  to  Ursell ’s  (1952)  solutions,  for  these 
or  related  topographies  and  if  so  how  may  they  be  found?  Second,  what  is 
the  appropriate  three  dimensional  edge  wave  solution  for  an  ’’almost”  linear 
topography?  Is  it  close  to  the  original  Stokes  solution  or  must  it  be 
close  to  the  high  frequency  limit  given  here? 

One  extension  of  this  work  would, be  towards  alongshore  variations  in 
topography.  If  the  topography  varied  periodically  along  the  cdastline 
(i.e.  In  y)  the  edge  waves  would  be  filtered  as  Floquet  waves, 
e.g.  Brillouin  (1953),  and  only  certain  frequencies  and  wavelengths  could 
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propagate.  Alternatively,  the  effect  of  a  small  irregularity  in  the 
otherwise  straight  coastline  could  be  studied,  e.g.  as  by  Fuller  and  Mysak 
(1977). 

The  extension  to  lower  frequencies  such  that  Coriolis  effects'  become 
significant  does  not  appear  to  be  appropriate  in  cdntext  of  a  three 
dimensional  theory  since  these  lower  frequencies  would  imply  shallow  water- 
long  wave  theory  to  be  a  good  approximation  in  the  real  ocean.  This,  then 
brings  up  the  question  of  the  applicability  of  two  dimensional  theories  to 
the  present  examples.  Case 'I-A  is  essentially  shallow  water  theory  but- 
the  remaining  asymptotic  solutions  do  have  some  depth  variation  as  dbes  the 
general  solution.  Thus  two  dimensional  theories  give  quite  different 
results  than  those  given  here  as  seen  in  appendix  B.  Even  case  fll-A, 
which  is  a  low  frequency  solution,  requires  very  small  h^  such  that  a-tyo 
dimensional  solution  nay  apply  for  x  <<  3^  eventually  (for  x  near  3^) 
vertical  variation  is  required. 

However  over  the  regions  of  interest,  l.e.  where  the  waves  have  not 
decayed  to  negligible  values,  even  the  asymptotic  solutions  do  not 
exhibit  appreciable  vertical  variations  for  r'tallstic  values  of  h*  and 

O 

Xjj  since  z  never  got  appreciably  different  from  1  (the  free  surface).' 

There  are  clearly  values  for  the  input  paraa.iters  where  these  vertical- 

variations  may  be  significant,  but  they  do  r.ot  appear  to  be  physically 
\ 

realistic,  e.g.  <  1. 
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APPENDIX  A 


The  following  computer  program  accept*  a  single  card  IJOB  as  the 
number  of  jobs  (input  data  seta)  and  then  several  "jobs"  each  with  h*.  h^j 
and  3^  (in  cm)  as  input.  The  card  IVEL  -  0  behind  IVEL  -  1  will  surpass 
velocity  output  while  the  reverse  will  provide  velocities  at 

"  0,0  and  “  0*0  (0.2)  1.0.  The  remaining  output  is  labeled; 

the  symbol  D  in  front  of  a  variable  means  the  dimensional  form  (in  cgs  units). 
The  symbol  HPX  is  h^  5^,  HOP  is  h^,  SIGMAO  is  aQt  etc..  Typical  running  times 
for  twenty  different  input  cards  is  about  2.5  seconds  on  a  CDC  CYBER  173. 
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DO  25  I  =  1*11 
X  =  T I -1 >  *XM2/ 1 y. 

IFd.EQ.lD  V  =  '.93***Mf 

ARG2  ~  CDSMIN2)  *ryp( -«*r  »*2*x /L?) 

h?(i  )  =  i  *  -  Accsm  a°s;  >"  2 

Y  s  X//‘-C 

PRINT  *  »A*‘  y=  ♦Y«CH  H2  =  «H2<IT 
I  F{  I  VEIL*  EO.i  )  GC  TO  2  5 
DO  255  *2  =  1*6 
DZ  =  H?(  I )  /  r». 

L  s  -T  17-1  )*^Z 
Xfjn  i  x/xm: 

XL2  =  X*L? 

I F  <  V  1.2  •"  T  •  lfO.)  FAC  =  a. 

I c  (XL2  *LE  •  rAC  =  •  rf  <-XL?  > /$IVH<  VC  > 

U(I*I7)  =  FAC*S  J'VMT  M2*  (1.  *2)  > 

YtI*iZ>  s  (K2/LD*FAC*SI5H<fv?*  (1  .♦?>  > 

\U«I7)  =  -Ci2/L2)*?»C*COSM(A,2*n  •♦Z>  > 
yRI  TE  (  6*  *•'  >  X*  XND*7  «,>LJ*lZ’>*Vtl*lZ>  *v»£  I*IZ> 

255  CONTINUE 
25  CONTIN'Uc 
29  CONTINI'T 

J  FT  »»py  » GT  •  0.5)  GC  Tt  55 
PRINT*.  5CM  C&5E  T  J.PX  *LT  .  '-.5 
XM-3  =  XR/H 
F  3  < 1 >  = 

DC  1M  0  =  2*11 
x  =  (.t-i  )*pi  n»  • 

F3TJ)  =  F  r  3 T  y ) 

C  PRINT  *  .Ajj  y=  *X  *6H  FT(J)=  *F3(J> 

101  CONTINUE 

DO  lit  J  =  111” 

TEST  =  CFMJ1  -P?X>*<  RTTO-»l)  -  HPX) 

ifttest  .lt.  c*j>  n  tj  112 

IFT  F 7  <  J  >  •  r-3X  )  5  =  J~1 
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-W«,s 


*h‘.>U?4.N  vftyyrLL 


7V17’  OPT=f  TRACE  WAN  TRAP 


FTN'  4.IM85 


:c<  FMj*  MPx  >  Gv  TO  58 

ii:  er.Tm:  ‘ 

1  H  ’ ,  t  T  *  » ';  M  Ci.T  "1*0  RrOT  .TO  r3  IN  THIS  RANGE 

GO  TO  <-<•,  , 

:  i:  V|  =  { J-1  )  *pj  /2  f.  ,  ‘ 

<h  -  .!♦"*  /;>'. 

V  *1?  ’  =  1*21  t 

**•*  =  UX  ♦  VL)/C. 

-  « »•:■■'  =  *V3(XM) 

IMvl  ’.v)  ">LT.  =  F3  (  1> 

!M  <L  .  T.  .,v)  FXL^  =  FF  3( XL ) 

y«=  *y«t  l.H  FF3<XM)  =  tFXM3  «  10H  FF3<XL>=  tFXl 
TL  -  -  n»X)  »  ( rXL7.  -  PPX) 

IF  (TL  .L  T.  )  XR  =  X '•■ 
l  r<  TL  .  r .  »•£  )  XL  =  VM 
Tet  AP$(  /P  -  XL)  .  L£»  ERKX-)  GO  TO  114 
13 

rRI,T»,A.r  rj  ^ALVJsC  ROUTINE  DID  NOT  PRODUCE  ROOT  1 
*c  =  tn»r>*>  xl=  «xl 

vO  "•*,  *• 

1*  •  =  (VH  ♦  XD/2. 

**’«  CO'  T!  Mi'  j  - 

K»I  r*t  *'M  »•;  fK  . 

IF  f  M.v  ..  GO  TC  *vB  :  .*  ,  »*  " 

ov  =  Si  «.T(y/Tfcf<(F>  >  -  i  *  V 

1. 3  s  V-  ' 

If  t ».  ?  .  L  T .  ••• )  '.0  TO  *  5 

*  ’  £  S£  H  T  ( I.  3  *  •  J  -  M ♦  ♦  2  >  ■ 

»  =  t>'»5H  u3=  »L  3  «5H  K3=  fK3rStt-$3s  «S3*6HV 
*•  *  A  =  -  T  /*.,  -  m  '  'i  • 


'I  = 

<k  = 

-i  _>  *  i  •  ’  r  ’ 

•*•*  =  (/,-i  ♦  X| 

-  «  V;T  =  -VT(V 


=  F  3  (  1> 


»FXH3  «  1LH  FF3  0XL)=  ,FXL3 


•  U 

r  v/v.'. 

“  c 
» 

,‘L7 

= 

-vy  i 

5  K  3  /  *1 

T!,’ 

-  i  f  .  p  I  C'i  A  C 

t  « I 

t  * « ^  -  :*••*  = 

T* 

=  T 5  /$  ? 

?(.r 

1>T’ 

=  r.  ♦"  uzf 

K  ; 

T  * »' -  nil  = 

>0 

~r  i  -  :  •  1 1 

v  s  f i-i »«x«3/i* . 

!F<  i  .{  r  .  J  *  )  y  =  :**93S*XM3 

=  (ASIN’{  EINC  *AMVA)*EXP(P.*«2*X/L3>>  -  GAMMA)  AM 

V  s  > />  ' 

f  /-  *Y  ir)H  i3=  ♦h2(I)  *  5 

If f  ? VtL. '  GO  TO  35 

or  ’5*  r  s  i,6  ■•'■ 

-  **?(]>/  r •  •  .  f 

2  =  -l?  ’-I  )*?Z  •  '- 

y*">  —  y / y  ' 

VL !  =  >*  1 3  .  ? 

FOL’.'T.  l'vi.)  FAC,  -  ?  *  , 

irrxL3.Lt  .  1  0.)  F4C  =  EX  P(  -  XL3)  /SIN {  M>  * 

I'UiIT)  =  FAC*S!N(v*tl,^Z>)  •  : 

V  l T  •  I  ?•)  =  C^/L3>*F4C*nN(  M*<1*»Z)> 

-  -<f./L3>  »f'AC*COS<  M*(l  •♦Z>)  * 

(t»oK>  x,x»o»z»u<itiz*),vn*iz>»w<iViz) 

’EG  Cu  <  7 1  f  'Uc 


PROGRAM  MAXWELL 


73/173 


f>|  Tr-'  M(  CE  ROSTRA* 


PTV  4.7*46r 


35  continue: 

GO  TO  ?5 

5  45  PRIt,T*»L3*  N*4UH  L3.LT.I' 

‘  55  CONTINUE 

290  ICOUNT  =  IC^UMT  *  1 

IF<ICPUNT  .LT.  1JO&)  TO  1 
’  END 


f 

r 

*  SYMBOLIC  REFERENCE  MAP  (R  =  !> 


Y  POINTS 


.1 

MAXWELL 

■A6LES 

;3 

AR61 

17 

OKI 

/C 

DK3 

53 

0L2 

56 

QM 

52 

0N2 

bl 

OS  2 

32 

DTI 

73 

QT3 

*7i 

CF  AC 

OOf 

ERRN1 

42 

FAC 

45 

FXL2 

20 

PXM1 

62 

FXM3 

64 

F  2 

>74 

G 

‘ns 

MM 

'10 

IIPX 

14  2 

H2 

n  r 

1 

'03 

I  JOB 

11 

'.63 

K1 

>67 

■K3 

>72 

L2 

>70. 

M 

,66 

N2 

746 

OLON2 

7"7 

SIG«AO 

759 

S2 

714 

'.TEST 

731 

»T1 

772. 

?3 

675 

W 

713 

X 

741 

XL1 

774 

-XL3 

71* 

'XMl 

TYPE  RELOC  *  T I 0\ 


REAL 

3  75  7 

A  A  7  2 

REAL 

3  754 

OK  J 

REAL 

T  7  2f, 

^Ll 

REAL 

2  7  A  7 

CL  3 

RrAL 

37  2?. 

0,1 

REAL 

3  73  ’ 

0S1 

REAL 

3-77) 

0S2 

REAL 

37C-S 

Dr2 

PEAL 

3736 

DE 

REAL 

37-* 

EKPN 

real 

3- 6 ’I 

V 

REAL 

2/21 

f  XL1 

REAL 

376’ 

FXL? 

REAL 

374  4 

rvvl2 

REAL 

1*1.337 

Fl 

REAL  ARRAY 

16411 

c  7 

REAL 

276- 

GAMS  A 

REAL 

2.7-4 

Hf*P 

REAL 

'77$ 

U1 

REAL  ARRAY 

h  2 ' 7 

u  1 

INTEGER 

3?:2 

I  COUNT 

INTEGER 

3-.7C 

IVF.L 

INTEGER 

3  712 

J 

REAL 

3  ?65 

KG 

REAL 

3  671 

LI 

REAL 

56  73 

L3 

REAL 

*664 

Ml 

REAL 

3723 

oLom 

REAL 

3n7E 

ri 

REAL 

3  724 

!»i 

REAL 

3764 

S3 

REAL 

3722 

TL 

REAL 

S7CC 

12 

REAL 

4454 

1) 

REAL  ARRAY 

153 16 

i 

REAL 

.57 15 

XL 

REAL 

376' 

yl: 

REAL 

27!  6 

t" 

REAL 

37  43 

y*u 

REAL 

REAL 

REAL 

RE  A  - 

Rr  A  L 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

P“AL 

INTEGER 

INTEGER 

INTEGER 

real 

REAL 

REAL 

REAL 

REAL 

REAL 

KrAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 


f  P  RA  Y«| 
array 


Af'RArl 
A  p  R  A  Y  ;| 


I 
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APPENDIX  B:  TWO  DIMENSIONAL  (SHALLOW  WATER)  EDOE  WAVES 

Consider  the  conditions  under  which  a  long  wavelength  free  surface 
gravity  wave  may  be  trapped  against  a  straight  coastline  by  a  topography 
which  varies  in  a  direction  normal  to  the  coastline.  Taking  x  in  this 
normal  direction  and  y  along  the  coastline,  the  linearized  shallow 
water-long  wavelength  equations  for  a  homogeneous  fluid  on  a  rotating 
earth  are  [c.g.  LeBlond  and  Mysak  (1978)]: 

3u/3t  +  g  3n/3x  ■  fv 

(1)  3v/3t  +  g  3n/3y  “  -fu 

3n/3t  +  H(3u/3x  +  3v/3y)  +  u3H/3x  »  0 

where  (u,v)  are  velocities,  q  is  the  free  surface  elevation,  H  is  the 
water  depth  and  f  is  the  Coriolis  parameter. 

Assume  solutions  which  propagate  along  the  coastline,  e.g. 

u(x,y,t)  •  U(x)  cos  (ky  -  ot) 

v(x,y,t)  -  V(x)  sin  (ky  -  ot) 

n(x,y,t)  -  E(x)  sin  (ky  -  ot) 

The  original  equations  then  reduce  to 

gE’  +  oU  -  fV  -  0 

(2)  -oV  +  gkE  +  fU  -  0 

-oE  +  H(U'  +  kV)  +  uH*  -  0 


C  2  ] 


Clearly  U  and  V  may  be  eliminated,  leaving  a  single  equation  on 
E,  e.g.  LeBlond  and  Mysak  (1978  -  p.  220) 

(3)  (HE*)'  +  |-q-2-~-f-2  -  k2H  -  ~  H’|e  -  0 

This  equation  does  not  force  a  particular  x  dependence  on  either  E  or 
H,  as  was  the  case  for  the  three  dimensional  equations,  since  here  E  arid 
H  are  coupled  while  previously  they  were  uncoupled.  Instead  there  are 
many  solutions.  If  for  example,  E  is  required  to  be  exponentially 
decaying  such  as  £Qexp  (**kc)  ,  equation  (3)  requires  H(x)  to 

be 


H(x) 


exp 


where  Hq  Is  an  arbitrary  constant  defined  In  general  either  by  choosing 
H(0)  »  0  or  U(0)  «  0  to  have  no  net  mass  flux  across  the .coastline, 
x  »  0.  If  H(0)  is  chosen  to  be  zero,  the  topography  is  given  by 


Certain  ranges  of  parameters  give  negative  topographies  which  must  be 
excluded.  Requiring  that  t  >  0,  but  recognizing  that  k  may  be  either 
positive  or  negative  yields  the  following  results: 

For  a  >  f,  k  >  l  gives  a  positive  concave  upwards  topography * 
exponentially  approaching  a  uniform  depth,  Hw  equal  to  '  2-f2)/g(k2-i^) 
-I  <  k  <  i  gives  a  positive,  convex  upwards,  exponentially  increasing 


[  3  ] 


depth,  -oX/f  <  k  <  -X  gives  a  positive,  concave  upwards 

topography  again  exponentially  asymptotic  to  and  finally  k  <  -o£/f 
yields  a  negative  topography. 

For  o  <  f,  k  >  X  and  -oX/f  <  k  <  X  both  yield 

negative  topographies,  but  -X  <  k  <  -oX/f  yields  a  positive, 

concave  upwards  topography  again  decaying  exponentially  to  the  asymptotic 
depth  H  while  k  <  -X  yields  a  positive,  convex  upwards 
exponentially  increasing  depth. 

The  choice  of  U(0)  equal  to  zero  requires  either  o  ■  -fk/X  , 

leading  to  a  constant  depth,  H  *  f2/gX2  ,  or  the  trivial  solution, 

E(0)  equal  to  Xero  and  thus  E(x)  equal  to  zero  everywhere.  This  choice  is 
then  rejected  for  this  particular  form  for  E(x) . 

In  general,  equation  (3)  may  be  considered  to  define  E  for  a  given 
H  or  alternatively  H  for  a  given  E.  While  the  first  view  is  more 
realistic  physically,  the  second  is  simpler  mathematically,  i.e. 


h-  -  &  e! 

H'  +  (e'  '  -  k2El  H  -  {——“I 

Q 

l  /  l  g  / 

may  be  solved  completely  for  H(x)  for  a  given  E. 


H(x)  •  Cq  exp 


1  [  £M(x) 


r-KlEfe}. 


E»(x)-  &  E(x) 


dx 


f  f2-o2 

l  8 


•1  exp[-  f  E"WfW-  d, 

'  I  ‘  E'W-^EO) 


h 


i  m 


E’(x)-  ~^  E(x) 


exp 


I  & 


-k2E(x) 


dx 


E'(x)-  f-  E(x) 


dx 


/ 


[  &  ] 


where  Cq  is  an  arbitrary  constant  to  be  chosen  to  have  zero  net  mass 
flux  across  x  ■  0. 

Although  this  format  may  seem  inappropriate  on  physical  grounds,  it 

represents  a  "transformation"  of  E(x)  into  H(x),  i.e.  trial  functions 

for  E  lead  directly  to  a  functional  form  for  H  and  in  a  sense  to  the  two 

dimensional  counterpart  to  the  three  dimensional  case  studied  in  the  main 

text.  The  alternative  view  of  specifying  H  and  solving  for  E  leads  to 

a  second  order  differential  equation  on  E  with  variable  coefficients  some 

of  which  vanish  when  H  is  zero,  i.e.  represent  singular  solution  points. 

A  number  of  solutions  in  this  format  have  been  generated,  e.g.  Ball  (1967) 

obtained  a  solution  in  terms  of  hypergeometrlc  polynomials  for  an 

exponential  depth  variation  of  the  form  H  (l-exp(-ax)) ,  Hidaka  (1976-b) 

o 

obtained  a  solution  in  terms  of  modified  Kathieu  functions  for  a  parabolic 

2  2  1/2 

depth  variation  of  the  form  Hq(1  +  (x  /a  ))  '  ,  Robinson  (1964)  used  a 

finite  width,  linear  shelf  terminating  at  a  discontinuity  to  a  constant 

depth  ocean  with  solutions  in  terms  of  Laguerre  functions,  Eckart  (1951), 

Mysak  (1968)  have  obtained  solutions  in  terms  of  Rummer  functions  or  Laguerre 

polynomials  for  a  linear  topography,  H(x)  •  yx,  and  Shaw  (1977)  has  obtained 

solutions  in  terms  of  (both)  Rummer  functions  for  the  case  of  a  linear 

topography  which  did  not  necessarily  break  the  surface,  H(x)  *  H  +  yx. 

o 

Other  solutions  could  be  generated  for  other  topographies  using,  the  standard 
functions  of  mathematical  physics,  e.g.  Bateman  (1954). 

In  a  sense,  however,  these  solutions  are  all  the  same  type; 
series  solutions,  using  the  method  of  Frobenius,  to  the  original  differential 
equations  with  the  singular  solution  (if  any)  at  the  origin  surpressed,  if 


necessary 
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